HW6 Interactive Plots with plotly()

1 Prepare a single data set

1. Life expectancy and income for the year 2015.

#Read data.
data0 <- read.csv ("/Users/qianshengliang/Dropbox/STA553/Week 6/LifeExpIncom.csv",
                   header=TRUE,sep=",", colClasses=c("NULL", NA, NA,NA,NA,NA,NA))
#DT::datatable(data0, fillContainer = FALSE, options = list(pageLength = 10))
#2015data
year.id=which(data0$year=="2015")
data2015 = data0[year.id,]
#range(data2015$pop) #11000 - 1400000000
#range(data2015$income) #623 - 120000
#range(data2015$lifeExp,na.rm=TRUE) #49.6 -83.8
#work data
data1 <- data2015
data1$pop2 <- (log(data2015$pop)-9)
#range(data1$pop2) #0.3-12
#write.csv (data2015, "/Users/qianshengliang/Dropbox/STA553/Week 6/data2015.csv")
DT::datatable(data2015, fillContainer = FALSE, options = list(pageLength = 10))

2 Make an interactive scatter plot

To display the association between life expectancy and income for the year 2015.

plot_ly(
  data = data1,
  x = ~income,  # Horizontal axis 
  y = ~lifeExp,   # Vertical axis 
  color = ~factor(country),  
  text = ~paste("<br>Country: ", country,
                "</br>Population: ", pop,
                "</br>Income: ", income,
                "</br>lifeExp: ", lifeExp),
  hoverinfo = "text",
  alpha  = 0.2,
  marker = list(size = ~(pop2), sizeref = 0.1, sizemode = 'area' ),
  type = "scatter",
  mode = "markers",
  ## graphic size
  width = 900,
  height = 700
) %>%
  layout(  
    ### Title 
    title =list(text = "Association between Income and Life Expectancy at 2015", 
                font = list(family = "Times New Roman",  # HTML font family  
                            size = 18,
                            color = "red")), 
    ### legend
    showlegend = FALSE,
    ## margin of the plot
    margin = list(
      b = 100,
      l = 100,
      t = 100,
      r = 50
    ),
    ## Background
    plot_bgcolor ='#f7f7f7', 
    ## Axes labels
    xaxis = list( 
      title=list(text = 'Income',
                 font = list(family = 'Arial')),
      zerolinecolor = 'black', 
      zerolinewidth = 2, 
      gridcolor = 'white'), 
    yaxis = list( 
      title=list(text = 'Life Exoetancy',
                 font = list(family = 'Arial')),
      zerolinecolor = 'purple', 
      zerolinewidth = 2, 
      gridcolor = 'white',
      range=c(40,90)),
    ## annotations
    annotations = list(  
      x = 0.8,   # between 0 and 1. 0 = left, 1 = right
      y = 42,   # between 0 and 1, 0 = bottom, 1 = top
      font = list(size = 12,
                  color = "darkred"),   
      text = "The point size is proportional to the population",   
      xref = "paper",  # "container" spans the entire `width` of the 
      #  lot. "paper" refers to the width of the 
      #  plotting area only. yref = "paper",  
      #  same as xref.
      xanchor = "center", #  horizontal alignment with respect to its x position
      yanchor = "bottom", #  similar to xanchor  
      showarrow = FALSE)
  )

The graph shows the relation between income and life expectancy of 193 countries/regions in 2015. The size of the points is proportioned to population and colors represent different contries.